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ABSTRACT 


This thesis presents a method for the determination of 
neutral gas densities in the lonosphere from rocket-borne 
measurments of UV atmospheric emissions. Computer models were 
used to calculate an initial guess for the neutral atmosphere. 
Using this neutral atmosphere, intensity profiles for the N, 
(0,5) Vegard-Kaplan band, the N, Lyman-Birge~-Hopfield band 
system, and the 0129724 line were calculated and compared with 
the March 1990 NPS MUSTANG data. The neutral atmospheric 
model was modified and the intensity profiles re-calculated 
until a fit with the data was obtained. The neutral 
atmosphere corresponding to the intensity profiles that fit 
the data was assumed to be the atmosphereic composition 
prevailing at the time of the observation. The ion densities 
were then calculated from the neutral atmosphere using a 
photochemical model. The electron density profile calculated 
by this model was compared with the electron density profile 
measured by the U.S. Air Force Geophysics Laboratory at a 


nearby site. 
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I. INTRODUCTION 


The earth's atmosphere has long been of interest to 
science. From explaining phenomenon such as the aurora 
borealis or northern lights to the study of weather and 
climate, knowledge of the atmosphere has many interesting and 
useful applications. In addition to the scientific community, 
the military considers an understanding of the atmosphere to 
be of extreme importance because of its profound effect on 
several critical systems. Specifically, the characteristics 
of the region of the upper atmosphere known as the ionosphere 
@re ©f particular significance. Among the features of 
interest, the electron density is of principal concern due to 
its role in high frequency (HF) radio communications, 
satellite communications and control, and over-the-horizon 
radar COTHR ). Another important aspect is the naturally 
occurring luminescence of the atmosphere called airglow. This 
1s significant because it offers a way to determine electron 
densities. Furthermore, airglowis the source of a background 
noise signal that affects infrared detection and tracking of 
ships, aircraft and ballistic missile launches by space 
systems. The electron density and airglow are not constant 
but vary due to several factors such as geographic location, 


time of day and solar activity. 


Currently, electron densities in the ionosphere are 
determined by actively sounding the atmosphere with radio 
waves. This is accomplished in the following way. A radio 
wave of a certain frequency is transmitted. This signal is 
reflected by the ionosphere when it encounters a unique 
electron density which corresponds to the transmitted 
frequency. This frequency in called the critical frequency. 
The reflected signal is detected by a receiver and the time it 
takes to make the round trip is recorded. By multiplying half 
this time by the velocity of waves in free space, one can 
determine what is known as the "virtual" height. This is not 
the actual height at which the reflection occurred. The 
calculated virtual height is actually greater than the true 
height of reflection. This is because the wave group velocity 
in the ionosphere is less than the free-space group velocity. 
Techniques exist for determining the true reflection height 
from the virtual height. By varying the frequency, an 
altitude-versus-density profile can be obtained. These 
soundings are conducted by instruments called ionosondes and 
while they are accurate, they have significant disadvantages. 
First, these are fixed sites with limited coverage and may not 
be located ina particular area of interest. Second, radio 
sounding is only capable of providing the lower half of the 
electron density profile from the ground. This is because 
although the same electron density may occur at different 


altitudes, the radio wave will always reflect off the first 


appropriate density it encounters i1.e., the one which is at 
the lower altitude. Finally, ionosondes do not provide a 
measurement of the airglow intensity or densities of other 
constituents of the ionosphere. 

As an alternative, one might suggest putting an ionosonde 
on a Satellite. This might correct the coverage problem but 
would then only provide the upper half of the density profile 
and would still not give information about the airglow or the 
additional densities of interest. Furthermore, due to its 
active nature (i.e., it transmits a signal), an ionosonde 
would necessarily be large with a large power requirement and 
thus not well suited for satellite deployment. 

An alternative solution is to place a spectrograph on a 
satellite which could record emissions from the airglow. 
These emissions could then be used to infer the densities of 
neutral and charged species from which the electron density 
can be determined. Because it 1S a passive instrument, a 
Spectrograph is less complex, requires less power and is 
smaller than an active system making it ideal for satellite 
deployment. One such instrument is the Naval Postgraduate 
School's middle ultraviolet spectrograph (MUSTANG). The 
MUSTANG was part of a scientific payload in two successful 
rocket-borne experiments: one conducted in March 1990 and the 
other in February 1992. This thesis analyzes the data 


collected from the 1990 flight. 


A. THESIS OBJECTIVES 

It is the purpose of this thesis to present a method of 
determining the densities of O, 0, and N, by comparing data 
with the model calculations. Additionally, the procedure for 
determining the densities of the charged species and electrons 


will be addressed. 


B. THESIS OUTLINE 

This thesis is divided into five chapters. Following the 
introduction, the second chapter gives background information 
necessary to understand the experiment and the analysis. 
Included in this chapter is a discussion of the earth's 
atmosphere, the airglow, and the atmospheric models. The 
third chapter is a brief description of the experiment 
including an explanation of the MUSTANG instrument's 
construction and how it works. Chapter IV contains the 
analysis of the data and the results. The final chapter 


discusses conclusions and suggestions for further research. 


ied bye BACKGROUND 
Before one can understand the analysis in Chapter IV, it 
is important to be familiar with the concepts and models on 
which this analysis is based. THis chapter provides the 
necessary background knowledge required to follow the 


analysis. 


A. THE EARTH'S ATMOSPHERE 

The earth's atmosphere extends from the surface out to 
approximately 1000 km. Figure 2-1 shows a temperature versus 
altitude profile of the atmosphere. One can see that the 
atmosphere 1s made up of several temperature gradients. These 
Gradients can be thought of as layers describing various 
temperature regions. The layer at the earth's surface is 
called the troposphere. Its temperature is influenced 
greatest by the surface temperature of the earth. This region 
1s characterized by a negative temperature gradient. The 
temperature continues to decrease until an altitude of 
approximately 10 km. At this point, called the tropopause, 
the temperature starts to increase with altitude. 

This next region is called the stratosphere. .The rise in 
temperature in the stratosphere 1s due to heating caused by 
the earth's ozone layer absorbing ultraviolet light. This 


region has received much attention in recent years because of 
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Figure 2-1 Temperature versus altitude profile of the earth's 
atmosphere. 


fe possible ozone-depleting effects of chloroflorocarbons 
(CFC's) in the atmosphere. These man-made chemicals, used as 
solvents and refrigerants, and in many industrial processes, 
are known to act as a catalyst in a reaction that destroys 
ozone. Without ozone, the earth would not be protected from 
the deadly effects of ultraviolet radiation and most life 
forms could not survive. The ozone heating peaks at the 
stratopause at an altitude of around 50 km. The next region 
is another negative temperature gradient layer called the 
mesosphere. Temperature decreases until the mesopause 1s 
reached at an altitude of about 86 km. Above the mesopause is 
the thermosphere where temperature increases until it becomes 
isothermal with increasing altitude. This isothermal 
Memperature is referred to as T, and varies between 1000 and 
1800K. This isothermal region is called the exosphere. The 
exosphere is characterized by extremely low densities and is 
the beginning of what one might call the "vacuum" of space. 
The region of interest for this thesis is the thermosphere. 
The major constituents of the thermosphere are molecular 
nitrogen, molecular oxygen and atomic oxygen. Their density 
profiles and the temperature profile in the thermosphere make- 
up what is referred to as the neutral atmosphere. Radiation 
from the Sun ionizes these species in a process called 
photoionization. The total ionization rate for a species at 


an altitude z, Q(z), is determined by the following equation: 


0 
Q*(z) =n*#(z) fot) 1, () e@n7t (A, 2) CH(®) yy (2m 
Ach 


where: 


n*(z) is the volume density of species a, 


o5(4) is the ionization cross section for species 4a, 


I,(») is the solar flux at the top of the atmosphere, 
*(z) is the optical depth at altitude z, 
CH(®) is the Chapman function. 


Aw, is the threshold wavelength 


The optical depth term reduces the total ionization rate at an 
altitude by accounting for the absorption in the atmospheric 
column above. The Chapman function accounts for the 
increasing column density as the sun moves from the zenith to 
the horizon. 

In addition to creating ions and free electrons, photons 


from the sun can create excited species and energetic 


electrons. These energetic electrons are called 
photoelectrons. Photoelectrons can cause additional 
lonizations and excitations. Also, molecules can dissociate 


as a result of photoabsorption or photoelectron impact. 
Chemical reactions also take place in the thermosphere which 
create other minor neutral and ion species. These species can 
also be ionized and excited. All of the above processes are 
called productions. Productions of one species result in the 
loss of others. Losses can also occur when dissociated 
molecules recombine or when excited particles relax back to 
lower energy states. The region formed by these reactions is 
called the ionosphere. 

The lonosphere 1S comprised of four regions D, E, Fl and 
F2. These regions are based on electron density. The density 
of each region changes with the time of day, season and solar 
activity. The D-region extends from an altitude of 50 to 90 
km, the E-region from 90-140 km and the F-region from 140- 
1000 km. The species in the ionosphere obey the one 
dimensional mass continuity equation (MCE) which is the 


subject of the following section. 


5. THE MASS CONTINUITY EQUATION 
The densities of all of the constituents in the 
1onosphere at a given altitude obey the following mass 


Sememmuity equation (MCE): 


On(Z) _ 2) . 
a =P-nL- (2-2) 


The term on the left-hand side of the equation is the change 
of density of a species at a given altitude with time. The 
first term on the right-hand side of the equation is the 
production processes discussed above. The second term on the 
right is the losses discussed above. The final term on the 
right is the flux term which accounts for diffuspomeor 
particles from outside the region. For most species, the 
entire MCE is required to calculate the density of a species 
at a given altitude. Some species are created and destroyed 
so rapidly that there is no time for them to diffuse to 
another altitude interval. As a result, the flux term can be 
ignored. In addition, for some species production and loss 
are in quasi-equilibrium, so the left-hand side of the 
equation is zero. Under these conditions, one can solve the 


MCE for the volume density: 


ayes a 
n(Z) ae (2-a 


This situation is called photochemical equilibrium (PCE). In 
the discussion below the atmospheric emission model uses PCE 
in the procedure that calculates volume emission rates. It 
calculates the volume emission rate from the volume excitation 
rate. This is done in the following way. If we assume that 
the entire density of an excited species relaxes giving off a 


photon, then the volume emission rate is the density times the 
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Binstein coefficient, where the Einstein coefficient Ais one 


over the lifetime of the species. The equation is: 


VE=nA (2=32)) 


Now substituting equation (2-3) for the density and setting L 
equal to A (because the losses are due to relaxation) we are 
left with the volume emission rate equal to the production. 
Since each species produced is assumed to relax and give off 
a photon, then the volume excitation rate is equal to the 
volume emission rate. There are some other losses such as 
quenching which must also be included but they do not change 
the above argument because they are simply dimensionless terms 
that reduce the magnitude of the production. Now that the 
atmosphere and the MCE have been discussed, another important 


topic, airglow, must be addressed. 


C. THE AIRGLOW 

As stated in the introduction, airglow is the naturally 
occurring luminescence in the atmosphere. Airglow is always 
present, day and night. imeeract  ~sehmemairgqiow for each of 
these time periods has its own name. Not surprisingly, they 
are the dayglow and the nightglow. 

One source of the airglow is the ionosphere. When the 
excited atoms and molecules of the ionosphere relax to lower 


energy states, they give off energy in the form of*a photon 


_— 


ee 


according to the equation: 3 


vas 


= yee 
E=hv ee ( 2 =e) 


The different species have several different energy levels 


Within a state. These slightly different energy differences 


result ina range of wavelengths. This range of wavelengths 
is called an emission band. When an intensity-versus- 
wavelength profile is made, the result is a spectrum. In a 


Simple case, only a single atomic or molecular emission is 
present for a given wavelength range. If this were the case, 
then the total intensity of the emission band can be 
determined by integrating over the wavelength range of the 
spectrum. Figure 2-2 shows a schematic example where one 
emission band produces the entire spectrum. Unfortunately, 
nature 1s much more complicated. A spectrum is the sum of the 
spectra of several processes. These processes produce 
emission bands that overlap or contaminate one another. 
Figure 2-3 shows the realistic case of a spectrum resulting 
from the sum of several contaminating emission bands. A 
technique called synthetic spectral fitting can be use to 
determine the intensity for the individual emission bands. 


This is described in Chapter III. 
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Figure 2-2 Simple schematic spectrum. 
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Figure 2-3 Complex schematic spectrum. 
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The 1990 MUSTANG data was collected during the day, 
therefore the emisSion bands used to analyze the data are 
dayglow emission bands. One of these emissions was the N, 
(0,5) Vegard-Kaplan (VK) band. The transition for the VK band 
system 1s 0 ume This is a forbidden transition because 
excitation from the ground state, xX, to the A state is 
forbidden by the selection rules of quantum mechanics. This 
does not mean that this transition cannot occur. It means 
that the A state cannot be populated by photoexcitation. The 
A state is populated in two ways: directly by photoelectron 
impact and indirectly by cascading from the C, B and W states 
(Meier, 1991). Another emission band observed in the data is 
N, Lyman-Birge-Hopfield (alll, - x22) : This is another 
forbidden transition. The a state is populated by 
photoelectron impact (Meier, 1991). Unlike the VK bands, this 
is the only production source for LBH. The final emission 
observed in the data was the 29724 line of atomic oxygen. 
This line is designated as 0129728. This is a ‘S - ‘P 
transition and it too is forbidden. 0('S) has different 
production and loss reactions for low-altitude (80-115 km) and 
high-altitude (above 115 km). In the low-altitude regime, the 
source of 0(*S) is the following reaction: 

Ono tM = 0, + M (2-5a) 
O. +M ~ 0, + O(‘S) (2-5b) 
where M is taken to be either N,, O0,, or O SO that the density 


[M] is equal to the sum of the densities (N.J, [0,], and [0]. 


i Eas) 


The sink for this altitude regime in addition to radiative 


relaxation is: 


O(7S) +02(a*A,) ~0+02. (2-6) 


The source reaction has a rate coefficient of 4.7x10°* (300/T)? 
cm’ s}, Because the proportion of 0, in the a Stateuimee 
unknown, the chemical loss term (Equation 2-6), is not 
included in the model calculations. This does not seriously 
affect the analysis because the low-altitude regime is below 
the region of interest. The source reaction also has little 
effect in the region of interest but was included for 


completeness. In the high-altitude regime the source 


reactions are: 


O+e*~0(1S) +e (2-7a) 
0; +e-0+0 (1S) We oy, 
O, +N-NO* +0 (1S) (2-Tep 


Equation (2-7b) has a rate coefficient of 2.1 x 10° em’ s? for 


a o(‘s) yield of 10% (B=0.1). Equation (2-7c) has a rate 


coefficient of 5 x 107°! cm s!. the sinks are the following: 
0(28) +0(1D) +hV..0, (2-8a) 
O(7S) -OCP) +h (2-8b) 


The Einstein coefficients for these two sinks are Assa7=1.06 gil 


and Ay977=0.045 gs), The above reactions and their 
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corresponding rate coefficients were taken from the review by 
Torr and Torr (1982). The next section describes the various 
models used in this thesis including how to use them, what 


they produce, and how they use the above information. 


D. THE MODELS 

This thesis uses three main computer models in its 
analysis. These models are the MSIS-83 model atmosphere 
(Hedin, 1983), MODATM, and the atmospheric emission model 
which was developed for this thesis. Each will be discussed 
below beginning with the MSIS-83 model. 

The MSIS-83 model, developed by Hedin (1983), computes a 
neutral atmosphere for a given set of conditions. When 
running this model, the user is prompted to enter several 
inputs. The first is IYD which is the year and date in the 
following format: year number and julian date i.e. YYDDD. 
Next is GLAT which is the latitude in degrees. Then the user 
1S prompted to enter GLONG which is the longitude in degrees 
from -180 to 180. After the longitude, one enters a number 
for F1O7A. This is the three month average of the F10.7 cm 
flux. The F10.7 cm flux is a radio signal produced by the sun 
which can be detected on the ground. It gives an indication 
of the sun's activity. Next, the user is prompted to enter 
Peeeewhich 1s the F10.7 flux for the day. The next. input is 


AP which is the A, index. This iS a measure of magnetic 
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activity. After AP, one 1s prompted to enter XLT which is the 
local time in hours on a 24 hour clock. The last input gs 
FRACXR which is the x-ray scaling factor and is another 
measure of solar activity. The values for F107, F107A, AP and 


FRACXR are available in NOAA publications of solar-terrestrial 


data. Once these parameters have been entered, MSIS-83 
computes the neutral atmosphere. The neutral atmosphere is 
produced in 2 km increments from 86 km to 450 km . The 


densities of O, 0», and N, as well as the temperature profile 
are stored in a file named YYDDDMSIS.DAT where the YY is the 
year and DDD is the julian date. The .DAT extension signifies 
that this is a data file. For this thesis, the file created 
was 9O0089MSIS.DAT. Any further reference to this file will be 
made as MSIS. To avoid possible confusion, the MSIS-83 model 
is referred to as MSIS3. This is shortened from the IDL 
version of MSIS-83 which is MSIS3.PRO. 

After the neutral atmosphere calculation is complete, the 
user is asked if he/she would like to calculate the 
photoelectron ionization rates for N,, O, and O, and the 
dissociation rate for N,. These calculations are performed by 
the Naval Research Laboratory (NRL) PEGFAC model (Strickland 
and Meier, 1982) which is called from MSIS3 but is a separate 
model. This creates a file 90089FLXM.PRT where the M 
Signifies that this file corresponds to the MSIS neutral 


atmosphere. Further references to 90089FLXM.PRT are made as 
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FLXM. If one wishes to modify the neutral atmosphere, one 
must use the next model called MODATM. 

The second model, MODATM, allows the user to modify the 
MSIS neutral atmosphere. The user runs MODATM after running 
MSIS. MODATM asks for the date in IYD format. The user is 
pee to change the following inputs: Tinf, T_120, T_86, 0_120, 
O2_120, N2_120 and the s parameter. The s parameter is used 
to calculate the temperature above 120 km. This parameter is 
not usually changed. Tinf is the isothermal temperature at 
the top of the temperature profile. T_120 and T_86 are the 
temperatures at an altitude of 120 km and 86 km respectively. 
The 0_120, 02_120, and N2_120 inputs are the densities for 
these species at an altitude of 120 km. When modifying the 
neutral atmosphere, the MSIS values are shown as a reminder to 
the user. For this thesis, only the 120 km densities and Tinf 
were changed. Changing T_86 and T_120 has much less affect on 
the profile than does changing Tinf. MODATM uses’ the 
densities at 120 km to re-calculate the density profiles. 
This re-calculation is affected by both the temperature at a 
gQiven altitude and Tinf. At lower altitudes, the effects of 
each are about equal. With increasing altitude, fTinf 
dominates. Therefore, only Tinf was adjusted for the analysis 
presented in this thesis. Once the changes have been made, 
MODATM calculates the new neutral atmosphere. It stores this 
information ina file called YYDDDNEU"V".DAT where the "V" is 


a version number to help distinguish between various runs. 


Jee, 


For example, the first MODATM neutral atmosphere file for this 
thesis 1s 90089NEU0.DAT. These files have the same format as 
the MSIS file. 

After the atmosphere 1s modified, PEGFAC is again used to 
calculate the photoelectron ionization rates and the N, 
dissociation rate. This 1s done while the MODATM is being 
tune This information is stored in a file called 
YYDDDFLX"V".PRT where the "V" is the same number as the "V" in 
the corresponding YYDDDNEU"V".DAT file. Further references to 
the MODATM neutral atmosphere and flux files will be made as 
NEU"V" and FLX"V". Both of these files have the same format 
as the MSIS3 files. Now that the calculation of a neutral 
atmosphere has been discussed, the procedure for calculating 
the intensity profiles for the VK, LBH, and 0129724 emission 
must be explained. 

The atmospheric emission model is a collection of IDL 
procedures that calculates the intensity profiles for the 
various emissions. The IDL codes are contained in the 
Appendix. The procedures that actually calculate the 
intensity profiles are SLNTINT.PRO and OISLNT.PRO. Although 
the profiles have been called intensity profiles, they are 
actually slant intensity profiles. A slant intensity 1s a 
column intensity at a given viewing angle. Unless specified, 
intensity will refer to a slant intensity. 

SLNTINT.PRO calculates the intensities for VK and LBH. 


One inputs the angle between the zenith and the instrument 
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viewing direction. This accounts for the difference in 
pathlength between the slant and vertical column intensities. 
For a detailed description of this see (Cleary, 1985). 
SLNTINT calls a procedure called GETVE.PRO. GETVE.PRO uses 
the following equation to calculate the volume emission rate 


pommetne (0,5) VK band: 


( (A3S+B3P+C3P/2+W3D) xB00x005) 
1+90x 10, (2-16) 


In equation (2-16), A3S, B3P, C3P, and W3D are the excitation 
rates for the A°*Z, Bll, C'JI, w3aA states, respectively; BOO is 
the branching ratio for the zeroth vibrational level of the A 
state; and Q05 is the Frank-Condon factor for the (0,5) 
transition. BOO takes into account that the A state is 
populated by direct excitation and cascade from the B,C, and 
W states. The second term in the denominator accounts for 
losses in the A state due to quenching by atomic oxygen. QO 
is the oxygen quenching rate and Ais the Einstein coefficient 
for the A state. Since the upper state for LBH is populated 
entirely by excitation to the a state, its volume emission 
rate 1S equal to the Al1P rate. The rates are passed into 
GETVE from the desired FLX"V" file by the procedure 
GETXRATE.PRO. Densities are pass into GETVE from MSIS or the 
NEU'V" file by GETDENS.PRO. SLNTINT then uses a procedure 


called VOLSLNT.PRO to convert the volume emission rates to 
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slant intensities. SLNTINT passes out an altitude array and 
the intensity arrays for VK and LBH. By plottimer the 
intensity arrays versus the altitude array, the intensity 
profiles for VK and LBH are produced. Now that the production 
of the VK and LBH intensity curves has been discussed, the 
production of the curve for O1I2972A must be addressed. 
OISLNT.PRO is the model that calculates the intensity for 
OI2972A. Prior to running OISLNT, one must run NOXIDIMGPRGe: 
This 1s a one-dimensional photochemical model which calculates 
the densities’ for 0, ' and N‘S, and the electron density 
(Cleary, 1985). The 0," and electron densities are stored in 
a file called YYDDDION"V".DAT and the N'S is stored in a file 
called YYDDDNOX"V".DAT. One inputs the viewing angle and a 
scaling factor into OISLNT. The procedure then asks the user 
to specify the desired neutral atmosphere (i.e., MSIS, NEUO, 
NEU1, etc.). Because PCE is assumed, OISLNT calculates the 
volume emission rate by dividing the production terms by the 
loss term and multiplying this by Ayg7)- The first prodWerroedm 
term comes from the reaction given by equations (2-5a) and (2- 
5b). The second production term comes from equation (2-7b). 
The third production term comes from equation (2-7c). The 
loss term is the sum of the two Einstein coefficients. ( To 
see the actual equations for the production terms see the code 
for OISLNT.PRO in the Appendix ). The procedure VOLSLNT 


converts this volume emission rate to a slant intensity. The 
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meeting Intensity array is plotted versus altitude giving 
the intensity profile for the OI12972A emission. 

Now that the background information has been covered, the 
experiment and analysis can be presented starting with the 


experiment in the next chapter. 
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lilies THE EXPERIMENT 


A. MUSTANG INSTRUMENT-CONSTRUCTION AND OPERATION 

The MUSTANG instrument has three distinct parts: the 
telescope, the spectrograph, and the electronics interface. 
Figure 3-1 shows a drawing of the MUSTANG instrument. A brief 
description of each section and their function follows. 

The first major component of the instrument is the 
telescope section. This is the collecting and focusing device 
for the instrument. The telescope is positioned off-axis from 
a 1/8th meter spherical mirror. Collimated light enters the 
telescope from the atmosphere and strikes the mirror which 
focuses the light at the entrance of the spectrograph. At 
this point, the light enters the mext section. 

The second section of the instrument is a modified Ebert- 
Fastie spectrograph. A 5mm x 140 wm entrance slit at the 
focus of the telescope mirror acts as the aperture of the 
spectrograph. The light then diverges from this point and 
falls on a 1/8th meter Ebert mirror. This mirror Wre= 
collimates the light onto the reflective diffraction grating. 
This 1200 lines/mm grating causes the light to reflect at an 


angle according to the following equation: 


asinO,=mA. CSEie 
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Schematic Diagram of MUSTANG Instrument. Major components 
are: (1) 1/8th m Telescope mirror, (2) Spectrograph enmance slit, 
(3) 1/8th m Ebert mirror, (4) Diffraction grating, (5) ITT Image 
intensifier, and (6) Hamamatsu Image detector. 


Figure 3-1 Schematic drawing of the MUSTANG instrument 
(electronic interface not shown) 


Zo 


The light then reflects off the Ebert mirror again and onto 
the ITT image intensifier in the 1800-3400A bandwidth. Due to 
the diffraction grating and the mirror, light of dijgi@ewen: 
wavelengths falls on the image intensifier in different 
positions. Once calibrated, it is the position that 
ultimately allows an intensity/wavelength relationship to be 
determined. The image intensifier is comprised of a quartz 
input window, cesSium telluride (CsTe) photocathode, two 
microchannel plates (MCP), a phosphor coated aluminum Screen, 
and a fiber optic o¥@tput “window, The image intensifier 
amplifies and converts the ultra-violet light to visible 
light. The photocathode generates a photoelectron when it is 
struck by a UV photon The electron is accelerated by an 
electric potential in the direction of the microchannel 
plates. When the electron strikes the MCP's, an electron 
avalanche occurs. The MCP acts as an electron multiplier 
producing approximately 15,000 electrons for every incoming 
electron. Another potential difference accelerates these 
electrons from the second microchannel plate to an aluminum 
screen. When the electrons strike the phosphor coating they 
cause the screen to give off:' visible light. The image 
intensifier is optically coupled to a Hamamatsu linear image 
sensor by its fiber optic output window and the image sensor's 
fiber optic input window. The image sensor uses 512 
photodiodes to convert the visible light intensity to a 


proportional voltage. 
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The third section of the MUSTANG is the electronic 
interface. This section uses an analog-to-digital converter 
amema first-in-first-out (FIFO) shift register to send the 
data to the telemetry section of the rocket. The data 
provided to the telemetry is sent to the ground via a RF 
Signal. For a detailed description of the interface section 


see (Quint, 1991). 


B. THE LAUNCH 

The MUSTANG instrument was launched on 30 March 1990 at 
1000 local time (1700 GMT) fromthe White Sands Missile Range, 
New Mexico. The launch vehicle was a NASA Terrier boosted 
Mmegeresralnt SOUNnding rocket. Figure 3-2 shows a basic diagram 
of the sounding rocket. The MUSTANG instrument was located in 
the experiment section of the payload. Other sections 
important to the experiment were the HV section which provided 
the power source and the T/M section which transmitted the 
data back to the ground. The payload reached an apogee of 
approximately 320 km. Data were recorded on both the up-leg 
and down-leg of the flight. This experiment resulted in the 
collection of airglow data in the form of 8000 spectra with a 
wavelength range from 1800A to 3400A. The spectra were taken 
between an altitude of 100 and 320 km. 

In order to determine the intensity of various emission 


bands, synthetic spectra of these emissions were fit to the 
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Figure 3-2 Launch vehicle block diagram. 
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spectra of the data. Synthetic spectra are spectra calculated 
for individual emission bands by models. The intensity of 
these spectra can be adjusted. By adding the synthetic 
spectra together and adjusting their intensities, fits for the 
data were achieved. Once these fits were obtained, the 
intensities for the various emissions were assumed to be the 
intensities of their corresponding synthetic spectra. 

The above process was performed by Walden (1991) using the 
N, (0,5) Vegard-Kaplan, N, Lyman-Birge-Hopfield, 299 ora omc 
oxygen and several other emissions. The 8000 spectra were 
averaged into 25 altitude bins. Each bin covered a 10 km 
altitude range. As a result, the intensity versus altitude 
profiles were determined for these emissions for the 1990 
flight. The profiles for the VK, LBH, and 01I2972A emissions 
were used in the atmospheric emission model to determine the 
neutral atmosphere. 

Now that the background information and the experiment 
have been discussed the data analysis and results can be 


addressed. 
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ive DATA ANALYSIS AND RESULTS 


A. ANALYSIS 
1. Theory 

The technique used to analyze the data is simple. One 
adjusts the inputs to the model until a satisfactory (2 
obtained. The neutral atmosphere that produces this fit is 
assumed to represent the conditions at the time of the 
observations. In this thesis, unless otherwise stated, the 
data refers to the intensity profiles from the analysis 
Denommed by Walden (1991). 

As a starting point, one uses the MSIS3 model to 
generate a neutral atmosphere and the atmospheric emission 
model to calculate an intensity profile for a "“Speq@uaae 
emission. The photoelectron excitation rates required by the 
model to produce intensity profiles for emission bands such as 
the Lyman-Birge-Hopfield, Vegard-Kaplan, and atomic oxygen 
2972K are created by the PEGFAC and stored in the FLXM data 
file. One then compares the intensity profile of the data 
with that of the model. If the two curves fit, then the 
atmospheric densities of O, O,, and N, (also written as [0], 
[0,], and [N,]) are the same as those calculated by MSIS3 and 


stored in the MSIS data file. 
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If the model calculation does not agree with the data, 
then one must modify the neutral atmosphere by adjusting Tinf 
ama/or [0], [0,], and [N,] at 120 km using MODATM. The model 
1S again used to generate an intensity profile and the effect 
that these changes have on the calculated profile is noted. 
re was found that an increase(decrease) in iUaloye 
increases(decreases) the scale height of the calculated 
emission while increasing(decreasing) (N, ] shifts the curve 
up(down) in altitude and increasing(decreasing) [0] shifts the 
Peorrle right(left). It was also found that changing [0, ] has 
little effect in the region of interest. These effects are 
consistent with the findings of Mack (1991). This process of 
adjusting the neutral atmosphere and comparing it with the 
data could lead to several sets of inputs which fit the data 
reasonably well in shape but not necessarily magnitude. 

Although the general shape of the profiles may agree, 
it may be necessary to multiply the curve by some constant to 
fine tune the fit. This is called a scaling factor. Scaling 
factors are generally small--usually less than ten. They are 
Significant because they have the effect of increasing the 
solar EUV flux used by MSIS3 and MODATM to predict the neutral 
atmosphere. For example, a scaling factor of 1.5 would 
increase the solar EUV flux by 50% which would have the end 
result of increasing the intensity by 50%. In addition, the 
Same scaling factor should be required to fit all other 


emission bands since the solar flux is the same for each. An 


Si 


inconsistency in the scaling factor indicates that the set of 
inputs used to produce the neutral atmosphere is suspect. The 
connection between intensity and the solar flux was described 
in Chapeer Tie 

Once a set of likely neutral atmospheres has been 
determined, then another emission is selected for examination. 
Intensity profiles for this emission are calculated using the 
atmospheric emission model and each of the different sets of 
inputs. The profiles for this second emission are compared 
with the data as before. Some of these calculated profiles 
will fit the data reasonably well while others will have the 
wrong shape, scale height, etc. and can be eliminated. One 
continues this process of selecting an emission, calculating 
the intensity curves using the remaining input sets and 
eliminating ones that produce curves that do not fit the data 
until only one combination remains. The density profiles 
stored in the NEU"V" file which corresponds to these initial 
conditions are those of the neutral species for that day and 
time. The density profiles of the various species of ions, 
excited species, and the electrons can then be calculated 
using NOX1DIM.PRO and the corresponding MSIS or NEU"V" file. 
These profiles are stored in 90089ION"V".DAT (90089NOX"V".DAT 
for the excited species) where again the "V" is either M for 


MSIS or the version number from MODATM. 
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Now that a general description of the theory has been 
discussed, the step-by-step procedure used to analyze the data 
can be presented. 

2. Procedure 

To begin the analysis, the MSIS3 program was run using 
the following values: AP=69, F10.7=186.5, ave Fl0.7=189.9 and 
X-ray scaling = 5 (Cleary, private communication, 1992). This 
Peeeemced a Tinf=1255K and the following densities at 120 km: 
[O]=8.8E10 em?, [0,]=4.4E10 cm?, and [N,]=3.2E11 cm’. It also 
created the data files MSIS and FLXM. Figure 4-1 is a 
comparison of the LBH curve generated by the atmospheric 
emission model using the above inputs with the data. As one 
can see, the calculated intensity profile for LBH shows a 
reasonable fit with the data both in scale height and in shape 
above approximately 120 km when scaled. The scaling factor 
necessary to get a good fit is 1.8. Although the model gives 
peemele information below 120 km, no attempt to fit the 
prediction to the data was made at lower altitudes for the 
following reason. Below, 20mm, tt 1S difficult to 
distinguish between the LBH bands and the brighter NO 
emissions. This caused the fit between the spectra and the 
data to be less certain. The intensity profile for the data 
obtained from the fit in the region below 120 kmis therefore 


less reliable. 
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Figure 4-1 Comparison of the data with the model intensity 
profile for LBH using MSIS values: .[0]_120=8. 8E10 lemme 
[0,]_120=4.4E10 cm’, [N,]_120=3.2E11 em”? and Tinf=1255K. 
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Evemetnough the above Erral produced a good fit, two 
other combinations of inputs were tried. These sets were 
suggested by Mack (1991). The purpose of this was to test 
whether either of these two neutral atmospheres would lead to 
a fit when compared with the data from Walden's analysis. To 
do this, MODATM was used to change the value of Tinf to 1500K 
and to reduce the MSIS3 [0] at 120 km by 27% and increase the 
MSIS3 [O,] and [N,] at 120 km by 5% and 30% respectively. This 
trial produced NEUO and FLX0O. These changes were made based 
on values suggested on page 82 of the report by Mack (1991). 
The effects of these modifications can be seen by the 
comparison of the LBH intensity curves of the atmospheric 
emission model with the data shown in Figure 4-2. The curve 
produced from these values has the wrong scale height and it 
also required scaling which was not needed by Mack. ; 

The other possible combination used Tinf=1700K and the 
input values for the neutral species densities at 120 km as 
suggested on page 83 of Mack (1991). These values represent 
a 48% reduction in [0], a 6% increase in [0, ] and the MSIS3 
value for [N, ] at 120 km. A run uSing these values was 
performed and the calculations were stored in NEU] and FLXl. 
Figure 4-3 shows the intensity profile generated by the 
atmospheric emission model feoreweueHh produced by this 
modification and the data profile. One can see that the curve 


produced using these inputs also does not agree with the data. 
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Figure 4-2 Comparison of the data with the model intensity 
profile for LBH using MODATM values: Tinf=1500K, -27% [0]_120, 
+5% [0,}_120 and +30% (N,J_120. 
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Figure 4-3 Comparison of the data with the model intens:t 
profile for LBH using MODATM values: Tinf=1700K, -48% Oa 20 
and +6% Kenmare ze : 
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When compared with the intensity profiles for the data from 
(Walden, 1991), neither set of inputs fit the data. The fact 
that this analysis rejected the sets of possible inputs 
suggested by Mack (1991) can be explained in the following 
way. The intensity profiles for the data were determined by 
fitting synthetic spectra to the data spectra. Mack (1991) 
used only part of the MUSTANG wavelength range while the later 
analysis by Walden (1991) used the full wavelength range of 
the instrument. Analysis over a smaller wavelength range, 
which includes fewer emission features, is much more sensitive 
to large fluctuations in the intensity than analysis over a 
larger range. One can imagine that if a fit was performed at 
the only large data spike in an otherwise even spectra, the 
fit would be poor for the whole spectral range. Because of 
this, the data profiles produced by Walden (1991) are 
considered more accurate. The input sets suggested by Mack 
(1991) do not fit the improved data curves. Therefore they 
were rejected. 

Since the original curve from MSIS/FLXM showed better 
agreement than either of the profiles from NEUO/FLXO and 
NEU1/FLX1, an attempt was made to fit the data using the 
Original Tinf=1255K while varying the neutral densities. 
Pursuit of a fit using NEUO/FLXO or NEUI1/FLX1 was not 
continued because the scale heights of these curves did not 


match the data. It was apparent that their Tinf values were 
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not correct and that the modification of the densities would 
Meeeresuit ina fit. 

Continuing the analysis, MODATM and the original MSIS3 
values for Tinf, [0,]_120, and [N,J_120 were used. Only half 
of the [0]_120 was used and NEU2/FLX2 was created. From 
Figure 4-4 which shows this LBH intensity profile and the 
data, one can see that this modification produced a curve with 
a peak intensity which was a little too low. As stated 
earlier, increasing [N,J_120 causes the peak intensity of the 
mode! calculation to occur at a higher altitude. The next 
trial makes use of this effect. To create the files NEU3 and 


FLX3, MODATM was run with all NEU2/FLX2 values except (N, ]_120 


was increased by 30%. The intensity profile produced with 
these values is shown in Figure 4-5. The effects of this 
change can be seen by comparing Figures 4-4 and 4-5. This 


last run, when scaled by 1.3, showed excellent agreement 
between the model and the data. Aw@eehis spommt, two input 
sets, MSIS and NEU3, still look capable of producing the 
correct neutral atmosphere so it was necessary to choose 
another emission and to continue the analysis. 

The next emission profile examined was the (0,5) 
Vegard-Kaplan band. When the model intensity profiles using 
MSIS and NEU3 were compared with the data, neither of the 
profiles for VK agreed well with the data as seen in Figures 
4-6 and 4-7. Both curves fit the peak intensity of the data 


well. However, moving above the peak, the agreement for both 
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Figure 4-4 Comparison of the data with the model intensity 
profile for LBH using MODATM values: Tinf=1255K and -50% [ Oma 
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Figure 4-5 Comparison of the data with the model intensity 
meoettle for LBH using MODATM values: Tinf=1255K, -50% [0]_120, 
and +30% [N,]_120. 
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Figure 4-6 Comparison of the data with the model intensity 
has gt for VK using MS{S values: Tinf=1255K, LO j— 120=8.8E10 
em~, [0.]_120=4.4E10 cm”, and [N,J_120=3.2E11 cm 
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Figure 4-7 Comparison of the data with the mode} intensity 
profile for VK using MODATM values: -50% [O]_120, +30% 


(N,]_120 amen tanr=1255K. 
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worsened due to a difference in the scale heights. 
Furthermore, in order to get the fit at the peak, it was 
necessary to scale the curve for MSIS by 4.8 and 2.7 for NEU3. 
Because the solar flux is the same for all emissions, the 
scaling factor should be the same regardless of the emission 
being examined. Therefore, any large variation in the scaling 
factor is an indicator that that set of inputs should be 
considered suspect. 

Because neither MSIS nor NEU3S lead to a good fijguee we 
order in which the emission features were fitted was changed. 
The intensity profile of the VK band is more reliable data 
than the profile for the LBH because the signal-to-noise ratio 
of the VK data was greater. As a result, reversing the order 
of the analysis to VK and then LBH is a better appaeach 
Furthermore, this ordering is superior because the model for 
VK is more complicated due to a greater number of production 
terms. This causes more variables to be constrained by 
analyzing the VK bands first. 

To begin the analysis using VK, a value of 3E-1l for 
the atomic oxygen quenching rate coefficient was determined by 
taking the average of the values reported by Piper et al. 
(1981), Thomas and Kaufman (1985) and De Souza et al. (1985). 
Also, an estimate of 0.125 for the zero vibrational level 
branching ratio was’ used. This was determined in the 
following manner. The first 8 vibrational levels (0 through 


Teor AS (no quench) from Figure 22 of (Cartwright , @iee 
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mremear to be of the same order of magnitude. Therefore, the 
branching ratio for the zeroth vibrational level was estimated 
to be 1/8 of the total relative number density of these states 
oo, 125. Using these values in the atmospheric emission 
model, previously rejected runs were re-evaluated. NEUO/FLXO 
and NEU2/FLX2 were again rejected because the curves showed 
poor agreement with the data (see Figures 4-8 and 4-9). On 
Mm@memmorner hand, the run using NEUL/FLX1 fit the data well 
except for a slight difference in the scale height which was 
only a little too high as seen in Figure 4-10. Although the 
calculations using MSIS/FLXM and NEU3/FLX3 were found tec be 
poor as stated above, they were the best two fits from the 
Original analysis attempt. Therefore, they were not 
immediately rejected. In addition, these runs, like the 
NEUI1/FLX1 trials, fit pretty well at the peak intensity but 
differed from the data elsewhere due to the scale height 
difference. 

One way in which the scale height can be increased is 
to increase the atomic oxygen quenching rate coefficient. 
Figure 4-1ll shows a VK intensity curve using 3E-11l1 and one 
uSing 9E-ll. This illustrates the effect that this change has 
Someeme VK intensity profile. This higher coefficient was 
suggested by Sharp (1971). Unfortunately, a change in the 
scale height is not the only effect that increasing the atomic 
oxygen quenching rate has on the VK profile. Figure 4-11 also 


shows that by increasing the rate, the peak intensity occurs 
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Figure 4-8 Comparison of the data with the mode! intensity 
profile for VK using MODATM values: Tinf=1500K, -27% [O}] zee 
+5% erp rose and +30% (N.]_120. 
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Figure 4-9 Comparison of the data with the model intensity 
profile for VK using MODATM values: Lines ok anda —50% 
woy 120. 
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Figure 4-10 Comparison of the data with the mode! intensity 
profile for VK using MODATM values: Tinf=1700K, -48% [O]_ ins 
and +6% [N,J]_120. 
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Figure 4-11 Effect of changing the atomic oxygen quenching 
rate coefficient on VK intensity profile. 
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at a higher altitude and the overall intensity drops. Tha 
drop in intensity iS most pronounced at the peak. The 
MSIS/FLXM, NEU1/FLX1, and NEU3/FLX3 runs were repeated using 
the higher quenching rate coefficient but the fits were again 
poor and they required an even larger scaling factor. 
Another way to increase the scale height of the VK 
emission profile is to increase Tinf as was discussed earlier. 
Taking this into consideration, one might guess that a 
slightly lower Tinf might improve the fit using NEU1/FLX1. 
Several attempts were made using Tinf=1600K, a quenching rate 
coefficient of 9E-11l and various combinations of neutral 
densities at 120 km. These trials were all variations of 
NEU4/FLX4 and none produced the desired effects (see Figures 
4-12 through 4-14). Scaling was attempted to fit the peak. 
Next, the quenching rate coefficient was changed back to 3E-]ll 
and better results were achieved as one can see from Figure 4- 
15. In this run, the peak intensity was too high so NEUS/FLX5 
was created using Tinf=1600k, half the MSIS3 [0]_120 value and 
increasing the MSIS3 [N,]_120 by 30%. This produced better 
scale height agreement but the altitude for the peak intensity 
was still a little too high as one can see from the intensity 
profile in Figure 4-16. The final run for VK used the 
NEU5/FLX5 values except the peak was lowered by lowering the 
[N,]_120 back to the original MSIS3 [N,]_120. This run was 


called NEU6/FLX6 and gave an excellent fit when scaled by 3.25 
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Figure 4-12 Comparison of the data with the model intensity 
profile for VK using MSIS values and changing the atomic 
eeeyoen quenching rate coefficient to q=9E-11l. 
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Figure 4-13 Comparison of the data with the model intensity 
profile for VK using MODATM values: Tinf=1600K, -50% [ O | eG 
and +30% (N-]_120; changing to G-=9RS J. 
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Figure 4-14 Comparison of the data with the model intensity 
profile for VK using MODATM value Tinf=1600 and changing to 
gq.=9E-11l ° 
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Figure 4-15 Comparison of the data with the mode! intensity 
profile for VK using MODATM value Tinf=l600K and changing to 


q-=3E-11 ° 
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Figure 4-16 Comparison of the data with the model intensity 
profile for VK using MODATM values: Tinf=1600K, -50% [0]J_120, 
+30% [N,J_120 and changing to q,=3E-1l1l. 


= 


(see Figure 4-17). This set of parameters was then used to 
fit the LBH data. This trial produced a good fit forutaeweee 
only at the peak intensity and required scaling the prediction 
by 1.5 (see Figure 4-18). A good fit over the full range of 
the profile could not be achieved. Due to the nature of the 
log/linear plot, any difference between the curve and the data 
will be greatest at the peak. Therefore, the best fit was 
determined to be the one that best matched the peak. 
Additionally, the signal is the largest at the peak so the 
spectral analysis that produced the data points is most 
reliable there. The difference between the scaling factor 
required for the VK fit and the LBH fit is addressed below. 

Previously, the analysis of VK was performed using an 
estimate for the branching ratio , denoted BOO. A more 
accurate value for BOO was determined by taking the ratio of 
the relative number density at the zero vibrational quantum 
number to the total relative number density for the N2 excited 
state A. The new value for BOO was found to be 0.25 and 
Figure 4-19 shows the effect that this change has on the VK 
intensity profile. This curve required a scaling Eae@temee: 
1.6 t6 daechievemassceocdmtams This is consistent with the 
scaling factor of 1.5 required for LBH curve as shown in 
Figure 4-18. 

The final emission used to analyze the data was the 


O012972A band. Using the model procedure OISLNT.PRO and 
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Figure 4-17 Comparison of the data with the model intensity 
profile for VK using MODATM values: Tinf=1600K, -50% nO} S120 
and changing to q,=3E-ll. 
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Figure 4-18 Comparison of the data with the model intensity 
profile for LBH using MODATM values: Tinf=1600K and -50% 
@ 4 MAG) 
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Figure 4-19 Comparison of the data with the model intensity 
profile for VK using MODATM values: Tinf=l1600K and -50% 
fo] 120: q-=3E-1l1 and changing to BOO=9.25. 
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NEU6/FLX6, a predicted intensity profile was produced for this 
emission and compared with the data. As discussed in Chapter 
II, the profile calculated by OISLNT.PRO is only accuratemro: 
altitudes above 115 skm.- Because of this, only the upper 
altitudes were considered when attempting to fit the model 
profile to the data. A fair fit was achieved when the curve 
was scaled by 2.1. This is shown by the comparison of the 
model O12972A intensity profile with the data in Figure 4-20. 
Even though the analysis for the 012972A emission produced 
only a fair fit, it 1s important to remember that this is the 
least reliable band. Because a problem with the telemetry 
caused some data to be lost in the region of this wavelength 
(Walden, 1991), there is a systematic uncertainty in the 29724 
intensity profile. With the systematic error and the error 
inherent in the manual fitting technique, the 2.1 Sc@aieae 
factor is acceptable. 

The above analysis indicated that the parameters from 
NEU6/FLX6 should correspond to the neutral atmosphere at the 
time of the launch. To test this, NOX1DIM.PRO (Cleary, 1986) 
was then run using these files and the profile for the 
electron density, [e ], was produced. This procedure uses the 
one dimensional MCE. A discussion of the results and a 
comparison of the calculated [e] profile with that of an 
ionosonde sounding performed by the Geophysics Laboratory (GL) 
at the time of the launch are presented in the following 


section. No comparison was made for the charged species 
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Figure 4-20 Comparison of the data with the model intensity 
Mofile for 012972A using MODATM values: Tinf=1600K and -50% 
foj_120. 


Sal 


because no independent density profiles of these species exist 
for this launch. 

In the course of writing this thesis, it was discovered 
that the W3D production rate was inadvertently omitted from 
the calculation for VK (see Equation 2-16). This omission 
caused a larger scaling factor to be required for all fits of 
VK. When the term was added, it was found that the best fit, 
NEU6/FLX6, fit without scaling. This is shown in Figure 4-21. 
Figures 4-22 and 4-23 shows two other figures for comparison 
to show that the ratio of the old scaling factor to the new 
scaling factor 1s approximately 1.6. This shows that the 
misSing term was incorporated in the scaling factor. This 
applies only to VK and an approximate 1.5 scaling factor is 
still required for LBH and = 3612972. This violates’ the 
requirement that the scaling factor for each emission be the 
same. If, for example, the scaling was accounting for a value 
for the solar flux which was too low, then the scaling factor 
would have to be the same for all emissions. The scaling 
factor is still important as an indicator of goodness of fit. 
Because the omission of the W3D term only resulted in a 
magnitude difference in all the VK plots, the above analysis 


would still have led to NEU6/FLX6 as the best fit. 
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B. RESULTS 

The above analysis lead to an estimate of Tinf=l1600K and 
a value for atomic oxygen density at 120 km of 4.4E10 cm”? 
which is half the MSIS3 calculated value. If NEU6 does in 
fact correspond to the neutral atmosphere at the time of the 
launch, then the [e ] curve calculated by NOX1DIM using NEU6 
would be expected to fit the [e] profile measured by GL. 
Figure 4-24 compares the NOX1DIM.PRO calculated [e] profile 
using NEU6 with the [e] profile from the GL ionosonde 
measurement (Bullett and Buchau, 1990). Figure 4-25 compares 
two NOX1DIM.PRO [e }] calculations with the ionosonde data: one 
uSing the MSIS neutral atmosphere and the other using the NEU6 
neutral atmosphere. This 1S intended to show whether any 
improvement can be seen in the calculated electron density 
profile by using the method describe above to determine the 
neutral atmosphere. As one can see, neither calculation 
agrees well with the GL data. The MSIS and the MODATM curves 
do seem to have the same general shape as the data from 90 to 
160 km and 90 to 280 km respectively. The conclusions drawn 


from the above results follow in the next chapter. 
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Figure 4-24 Comparison of the ionosonde data with the NOX1DIM 


calculated [e] profile using MODATM values: Tinf=1600K and 
moos [O}] 120. 
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Figure 4-25 Comparison of the ionosonde data with the NOX1DIM 


is [e ] profiles using MSIS values and using NEU6/FLX6 
values. 
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V. CONCLUSIONS 
A. SUMMARY 

It was the aim of this thesis to develop a new technique 
for determining the densities of OO, 0, and N, in the 
lonosphere. These species make-up the neutral atmosphere. It 
was hoped that by comparing calculated intensity profiles with 
the data, one could adjust the inputs until the model 
intensity profiles fit the data. The neutral atmosphere 
calculated by the model would then correspond to the neutral 
atmosphere at the time that the data was taken. This could 
then be used to determine the densities of the charged species 
and electrons. 

The analysis method presented in this thesis did not, at 
first glance, appear to provide a viable approach for 
determining the neutral atmosphere or the densities of the 
charged species and electrons. This is somewhat misleading 
for several reasons. First, the data from the GL ionosonde 
measurement was taken 310 km north of the White Sands launch 
Site. Normally, this would provide an accurate electron 
density profile for White Sands but at the time of the launch 
there was higher than normal geomagnetic activity. This 
resulted ina Traveling Ionospheric Disturbance or TID which 
distorted the ionosphere causing lower peak densities to occur 


at higher altitudes. These effects were reported by Bullett 
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and Buchau (1990) but could not be taken into account in the 
models used in the analysis. Also, since the ionosonde 
measurement was not taken in White Sands, the effects of the 
TID at the launch site could only be estimated and thus the 
fe ] is also an estimate. Second, the experiment lost some 
data in the 2972A wavelength region due to a problem with the 
telemetry circuit (Walden, 1991). Third, only the three 
emission bands were used. Using other prominent emissions 
would have helped the analysis. One emission that would have 
been particularly useful is the Second Positive band of 
molecular nitrogen. An analysis of this prominent band could 
not be performed because a synthetic spectrum model for the 
data currently does not exist. Another method for determining 
the intensity profile for the Second Positive suggested by 
Walden (1991) proved to be unworkable. Had the O12972A 
emission been more reliable and if the Second Positive band 
had been available, a better calculation might have been 
possible uSing the above analysis method. Finalise the 
scaling factor of approximately 1.6 which was necessary to 
achieve a good fit for LBH and 0129724 is not understood. It 
suggests that there is uncertainty in the models for the LBH 
and OI12972A emissions. The excellent fit of VK also suggests 
that the data for both LBH and 012972A is less reliable than 
Originally believed. One last conclusion that can be drawn 
from this thesis is that the excellent fit for VK and LBH 


(once scaled) in scale height indicates a Tinf=1600K which is 
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higher than the Tinf calculated by the MSIS3 model. It also 
is different than the 1500K and 1700K suggested by Mack 
(1991). The Tinf=1600K resulted in much better scale height 
agreement for the intensity curves. 

The method described in this thesis did not produce a 
neutral atmosphere which led to an electron density profile 
that fit the Geophysics Laboratory data. Although this is 
true, this approach cannot be dismissed. Improving the 
problem areas discussed above should allow this method to 
successfully calculate the neutral atmosphere as well as the 


charged species and electron density curves. 


B. SUGGESTIONS FOR FOLLOW-ON RESEARCH 

As pointed out in the summary section, there are several 
areas in which improvements can be made to better test the 
viability of this method for analyzing the ionosphere. First, 
improvements oes made for the 1992 launch which should have 
increased the accuracy of the data collected. One such 
improvement was better telemetry electronics which eliminated 
the data drop problems that hurt the information collection in 
1990. Using this better data, the synthetic spectrum fitting 
performed by Walden (1991) should be repeated to get the 
intensity profiles for VK, LBH and OI2972A. The elimination 
of the drops in the data should result in a more accurate 
intensity profile for OI2972A. Additionally, the model for 


OI12972K should be extended to include the lower altitude 


7: 


profile of this emission. Second, a synthetic spectrum should 
be developed for the Second Positive band. As was suggested 
by Walden (1991), the known contributions of N, could be 
subtracted from the data and the remainder could be assumed to 
be the Second Positive spectrum. This should now be practical 
with the better 1992 launch data. Third, other emissions 
should also be used in the fit process. Using several 
emissions should lead to better results. Fourth, the scaling 
requirement should be more fully investigated. Possible 
refinements in the models for LBH and 012972 as well as the 
better data for the 1992 launch should be able to eliminate 
these. Lastly, an effort should be made to fully integrate 
the atmospheric emission model, the MSIS3 model, MODATM, and 
NOX1LDIM into a single program that uses numerical methods to 
automatically determine the fit. This would be necessary to 
make this method practical for analyzing the large amount of 
data received from a satellite-based spectrograph. If the 
above changes are implemented, this method and the data from 
satellites could be used to make timely ionospheric profile 


calculations possible for any point on the globe. 
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APPENDIX 
The Appendix contains the IDL code for the Atmospheric 


Emission Model. 


*;THIS PROCEDURE GETS THE 0,02,N2 
‘DENSITIES 

*;AND CORRESPONDING ALT FROM EITHER 
;MSIS.DAT OR NED’V’ .DAT 

‘FILES FOR USE IN GETVE.PRO. 
*;WRITTEN BY KENNETH L. MCELROY 

715 MAY 1992 


PRO GETDENS, OALT, ODEN, O2DEN, N2DEN, TEMP . 
BILE=’ ' 

READ, ’WHAT DATA FILE DO YOU WANT TO USE? (I.E. 9O0089MSIS.DAT OR 90089NEU"V".DAT) ’ 
,EILE 

OPENR, 2, FILE 

B=ASSOC (2, FLTARR (183) ) 

OALT=B (1) 

TEMP =B (2) 

ODEN=B (3) 

IZ2DEN=B (4) 

N2DEN#=B (5) 

CLOSE, 2 

RETURN 

END 
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;THIS PROG READS THE ALT AND VOL 
*;PROD RATES A1P,A3S,B3P,C3P, AND W3D 
;FOR N2 FROM FLX’V’.PRT OR 
*;FLXM.PRT FILES 
*WRITTEN BY K.L. MCELROY 
723 APR 1992 
*;MODIFIED 06 JUNE 92 
PRO GETXRATE,ALT,A1P,A3S,B3P,C3P,W3D,0I1 

file=’ ’ 

read, ’What file do you want to use? (FLX"V".PRT OR FLXM.PRT)’, file 

OPENR,1, FILE 


;THIS SECTION GETS ALT AND A3S 
:;FOR N2 
STRI1=’ ? 
WHILE STRPOS(STR1,’VOLUME PRODUCTION RATES’) EQ -1 DO READF,1,STR1 
FOR J#1,5 DO READF,1,STRI1 
ALT#FLTARR (51) 
A3S#FLTARR (51) 
Rl#FLTARR (11) 
FOR J=#0,50 DO BEGIN 
READF,1,R1l 
ALT (J) =R1 (0) 
A3S (J) =R1 (10) 
ENDFOR 


*THIS SECTION GETS A1P,/ Bapeeee 
; AND W3D 
STR2=’ ¢ 
FOR J#1,7 DO READF,1,STR2 
AlP=FLTARR(51) 
B3P=FLTARR(51) 
C3P=#FLTARR (51) 
W3D#FLTARR (51) 
R2@FLTARR (5) 
FOR J#0,50 DO BEGIN 
READF,1,R2 
AlP (J) #R2 (4) 
B3P (J) #R2 (1) 
C3P (J) #R2 (2) 
W3D (J) =R2 (3) 
ENDFOR 
WHILE STRPOS(STR2,’OI’) EQ -1 DO READF,1,STR2 
FOR J=1,3 DO READF,1,STR2 
OI#FLTARR (51) 
FOR J=#0,50 DO BEGIN 
READF,1,R2 
OI (J) =R2 (3) 
ENDFOR 
CLOSE,1 
RETURN 
END 
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;THIS PROCEDURE USES GETXRATE.PRO AND 
;GETODEN.PRO TO GET VOLUME EXCITATION RATES 
iFOR VK, LBH, AND SECPOS FOR SLNTINT.PRO 
fWRITTEN BY K.E. MCELROY 
715 MAY 1992 
;MODIFIED:14 SEP 92 

PRO GETVE, ALT, VVK, VLBH, VSECPOS, DENOM 

GETXRATE,ALT,A1P,A3S,B3P,C3P,W3D 

GETDENS, OALT, ODEN 


=.5 pfH1S 1S L/LIFETIME OF A STATE(S*“=-1) 
BOO=.125 ;THIS THE BRANCHING RATIO 

Q0=3E-11 sTHIS IS © QUENCHING COEFF. (CM*3 S%*-1) 
O05=1.69E-1 ;THIS IS THE FRANK-CONDON FACTOR FOR 0-5 
O=INTERPOL (ODEN, OALT, ALT) 
VLBH=A1P 


VSECPOS=C3P/2. 

DENOM= (1 .+Q0*0O/A) 

VVK= ( (A3S+B3P+C3P/2.)*B00*Q05) /DENOM 
RETURN 

END 

@GETXRATE 

@GETDENS 


(he 


;THIS PROCEDURE CONVERTS VOL EX 

*;FROM GETVE.PRO TO SLANT INTEN 
;FOR VK, LBH, AND SECPOS OF N2 

*;WRITTEN BY K.L.MCELROY 

723 APR 1992 

;MODIFIED: 15 MAY 92 


pro slntint,att0,alt,svk,slbh, ssecpos, denom 
getve,alt,vvk,vlbh, vsecpos,denom 

att=alt*O0+att0 

delalt=alt (0) -alt (5) 

tophvk=abs (delalt/alog(vvk (5) /vvk(0))) 
topvevk=vvk (0) *exp(-2./tophvk) 

tophlbh=abs (delalt/alog(vlbh (5) /vlbh (0) )) 
topvelbhe=vlbh (0) *exp(-2./tophlbh) 

tophc3pisabs (delalt/alog(vsecpos (5) /vsecpos (0))) 
topvec3pi=#vsecpos (0) *exp(-2./tophc3pi) 

volsint, vvk,alt,att,topvevk,tophvk,svk 
volsint,vlbh,alt,att,topvelbh,tophlbdh,slbh 
volslnt,vsecpos,alt,att,topvec3pi,tophc3pi,ssecpos 
svk=svk*le-6 

slbh=slbh*le-6 

ssecpos=ssecpos*le-6 

return 

end 

Rgetve 
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PR OISLNT, ATTO,ALT183,S0OI1 
ATT=FLTARR (183) 
ATT#ATT*0O+ATTO 

A2972=0.045 

A5577=1.06 

B=.1 

K2=2.1E-8 

K3=5Eb-11 

Va! f 


READ, ’WHAT VERSION OF FILES DO YOU WANT?’,V 


OPENR,1,/’ 90089FLX’ +V+’ .PRT’ 
STR1=’ 


*THIS PROCEDURE CALCULATES THE SLANT 
s INTENSITIES FOR 0O1(2972). 

*WRITTEN BY K.L. MCELROY 

sll JULY 1992 

*MODIFIED: 16 SEP 92 


*THESE ARE THE EINSTEIN COEFFIENTS 


*THIS IS THE BRANCHING RATIO 


*THIS SECTION GETS THE FIRST 
PRODUCTION FACTOR Pl 


WHILE STRPOS(STR1,’VOLUME PRODUCTION RATE’) EQ -1 DO READF,1,STRI 
OI’) EQ -1 DO READF,1,STR1 


WHILE STRPOS(STR1,’SPECIES 
FOR J#1,3 DO READF,1,STR1 
ALT=FLTARR (51) 
O1S=FLTARR (51) 
R1=FLTARR (11) 
FOR J#0,50 DO BEGIN 
READF,1,R1 
ALT (J) =R1 (0) 
01S (J) =R1 (3) 
ENDFOR 
CLOSE,1 


OPENR,1,’ 90089ION’ +V+’ .DAT’ 
A=ASSOC (1,FLTARR (183) ) 
QO2P=FLTARR(183) 

ELEC#FLTARR (183) 
%2P=REVERSE (A(5) ) 
ELEC=REVERSE (A (6) ) 

CLOSE,1 : 


OPENR,1, ’90089NOX’ +V+!’ .DAT’ 
K=ASSOC (1,FLTARR (183) ) 
N4S=FLTARR (183) 

N4S=REVERSE (A(1)) 

CLOSE, 1 
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*THIS SECTION GETS THE 
*DENSITIES FOR O2+ & ELEC 


*;MAKES O2P AND ELEC ARRAYS GO 
*;FROM HIGH TO LOW 


THIS SECTION GETS N4S DENSITY 


*;MAKES ARRAY GO FROM HIGH TO LOW 


*THIS CONVERTS O1S FROM 
7A FLTARR(51) TO FLTARR (183) 


ALT183=FINDGEN (183) *2.+86 

ALT183=REVERSE (ALT183) 

O1SDEN=INTERPOL(01S,ALT,ALT183) 
*;THIS SECTION GETS 0,02 & N2 
*;DENSITIES AND TEMP NEEDED FOR 


;P4 
GETDENS, OALT, ODEN, O2DEN, N2DEN, TEMP 
ODEN#=REVERSE (ODEN) ;TBESE COMMANDS MAKE THE 
O2DEN#=REVERSE (O2DEN) ;DENSITIES GO FROM HIGH TO LOW 
N2DEN=REVERSE (N2DEN) 
M=ODEN+02DEN+N2DEN 
T=REVERSE (TEMP) ;MAKES TEMP GO FROM HIGH TO LO 


K4=4.7E-33* (300/T) *2 
*THIS SECTION CALCULATES THE V 
;FOR 012972 
P1*01SDEN 
P2=K2*B*O2P*ELEC 
P3#K3*O2P*N4S 
P4=K4*ODEN“%2*M 
L=A2972+A5577 
VOI=(P1+P2+P3+P4) /L*a2972 
;THIS SECTION CALCULATES THE 
;SLANT INTENSITY FOR 012972 
DELALT#ALT183 (0) -ALT183 (5) 
TOPHOI=ABS (DELALT/ALOG (VOI (5) /VOI (0) )) 
TOPVEOI=VOI (0) *EXP (-2./TOPHOT) 
VOLSLNT, VOI, ALT183,ATT, TOPVEOI, TOPHOI, SOI 
;THIS CONVERTS TO RAYLEIG&S 
SOI=SOI*1E-6 
RETURN 
END 
@GETDENS 
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